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Abstract 

The dynamics of a tracer particle in a glassy matrix of obstacles displays slow complex transport as the free volume 
approaches a critical value and the void space falls apart. We investigate the emerging subdiffusive motion of the test 
particle by extensive molecular dynamics simulations and characterize the spatio-temporal transport in terms of two- 
time correlation functions, including the time-dependent diffusion coefficient as well as the wavenumber-dependent 
intermediate scattering function. We rationaUze our findings within the framework of critical phenomena and compare 
our data to a dynamic scaling theory. 
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1. Introduction 

The structural dynamics in glass-forming liquids slows 
down by orders of magnitude upon cooling or compres- 
sion eventually leading to a quasi-arrested state or a dis- 
ordered solid. Conventionally the dynamics close to 
the glass transition slows down uniformly as manifested 
in the divergence of a single time scale characterizing 
the slowest process in the system, referred to as the a- 
process. A coherent theoretical picture for a variety of 
phenomena has emerged by the mode-coupling theory 
(MCT) of the glass transition fTH, developed by Wolf- 
gang Gotze and collaborators in the last 25 years. In 
particular, this approach has explained how a two-step 
relaxation process with a nontrivial two-time fractal re- 
sults from the equations of motion close to a bifurca- 
tion point. The mechanism encoded in these equations 
to yield power laws with exponents that are not simple 
fractions has no precedence in other areas of physics and 
yet has been identified as generic in nonUnear integro- 
differential equations. 

The mathematical properties of the MCT solutions 
are well understood for single component liquids ^ 
and mixtures |3|; in particular, it has been shown that 
the solutions allow a representation in terms of a contin- 
uous distribution of relaxation rates Q. The success of 



the mode-coupling theory is based on the fact that it not 
only explains how the slow dynamics emerges in prin- 
ciple, but that it also has provided a series of testable 
predictions which challenge the physical intuition ob- 
tained so far. To mention just a few, MCT coherently 
explains the physics of suspensions of hard-sphere col- 
loids |l5]|6l|7l, the reentrant phenomenon in 'attractive 
colloidal glasses' p8','9l, a pronounced minimum in the 
light-scattering spectra of supercooled liquids ITOl [TTl 
[T2l . or the composition dependence of the structural 
relaxation in mixtures |TT, TT]. For many years only 
the standard mode coupling equations have been dis- 
cussed and attempts to go beyond have been focused on 
extension for non-spherical molecules ifTSl [T6l [TTl [TSl 
and more recently to include shear in colloidal suspen- 
sions |fT9ll20l . 

The overall success of MCT is encouraging to inves- 
tigate more complex glass-transition scenarios to obtain 
a deeper insight into the nature and quality of the ap- 
proximations involved. A particularly interesting candi- 
date are strongly size-disparate mixtures where struc- 
tural arrest may occur in several steps. There, the 
majority component of large particles can undergo a 
glass transition characterized by a frozen disordered 
structure, whereas the minority species of smaller size 
meanders through the emerging network of channels. 
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Neutron scattering experiments and molecular dynam- 
ics simulations indicate a significant separation of time 
scales for sodium silicate melts also at finite concentra- 
tion ETII . which could be rationalized within standard 
MCT calculations (2T\ later. A similar splitting of relax- 
ation times has been found also for size-disparate soft 
spheres ll23l and Yukawa mixtures ll24ll . The dynamics 
within a frozen matrix or nanoporous medium has been 
studied only recently |,25. ,26.1 corroborating an intrigu- 
ing interplay of several dynamic transitions predicted 
within an extension of the mode-coupling theory for the 
dynamics within a disordered matrix l27il28L The stan- 
dard MCT of mixtures appears to qualitatively describe 
many aspects of the splitting of the dynamics (29), yet 
ultimately predicts that the structural arrest of the ma- 
jority and minority particles occurs at the same criti- 
cal point. The observed peculiarities are rationalized 
as a precursor phenomenon for the two-step freezing 
not contained in the standard MCT. Krakoviack's exten- 
sions to disordered matrices treats the frozen obstacles 
and the fluid on unequal footing 127.. 30] thereby allow- 
ing for multiple transitions. Although both approaches 
appear to give a satisfactory qualitative picture, theo- 
retical issues remain that are poorly understood. First, 
the freezing in the disordered solid is accompanied by 
a divergent length scale suggesting that the transition 
is driven by an entire hierarchy of length scales rather 
than the Lindemann length for caging. It appears that 
MCT in its current form is ill-suited to deal with this 
phenomenon 1 1 1. This manifests itself in 'spurious long- 
time anomalies' |30| that, if treated correctly, lead to a 
change in the exponent of the anomalous transport ob- 
served in the mean-square displacement |31 1. 

An alternative approach dealing with the dynamics of 
a single tracer in a disordered matrix is within the frame- 
work of critical phenomena. The localization transition 
of the tracer is due to an underlying percolation transi- 
tion of the void space as the density of the matrix is al- 
tered. As is well known, this geometric problem leads to 
a self-similar distribution of clusters Il32l entailing a se- 
ries of scaling predictions. It has been suggested that the 
dynamics of the Lorentz model shares the same univer- 
sality class with a random resistor network with power- 
law distributed conductances f33l [3411351 . The critical 
dynamics of the Lorentz model in the vicinity of the 
percolation threshold can be described within a scaling 
ansatz for the van Hove correlation function [36 1. Fur- 
thermore, corrections to scaling due to irrelevant scaling 
variables are not negligible and have been included re- 
cently 137] Ea. 

In this paper we investigate the Lorentz model by ex- 
tensive computer simulations and show that the mean- 




Figure 1 : Typical configuration for a random array of obstacles where 
the void space for a tracer is close to the percolation threshold. The 
tracer particles (highlighted) displayed have equal radius as the obsta- 
cles and are confined to the void space. 

square displacement becomes anomalous with an ex- 
ponent differing from a simple fraction. The vanish- 
ing of the diffusion coefficient upon approaching the 
localization transition is connected to the subdiffusive 
motion and the divergence of the correlation length or 
mean cluster size. We discuss, in particular the mo- 
tion of tracers that are confined to the percolating clus- 
ter in terms of the mean-square displacement, as well as 
the time-dependent diffusion coefficient averaged over 
all initial conditions of the tracers. We then elucidate 
the transport as response to an alternating external field 
in terms of the frequency-dependent conductivity and 
susceptibility. Furthermore, we characterize the space- 
resolved dynamics in terms of the intermediate scatter- 
ing function, which is accessible to scattering experi- 
ments in principle. Its long-time limit, known as the 
Lamb-MoBbauer factor or nonergodicity factor, reveals 
the trapping of particles inside the finite clusters. 

2. The Lorentz model 

A simple model for transport in disordered solids cap- 
turing all the relevant ingredients for complex trans- 
port was introduced by Lorentz already in 1905 f39\. 
There a single tracer or equivalently a system of parti- 
cles which do not interact among themselves traverses 
a course of immobilized obstacles. The position of the 
tracer is excluded from the space occupied by the obsta- 
cles thus confined to the void space. In the simplest case 
hard spherical obstacles are distributed independently 
and randomly in space and consequently the number of 
obstacles per volume V characterizes the emerging 
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structure completely. Then a dimensionless control pa- 
rameter n* - Ncr'' IV can be introduced, where cr is the 
distance of the hard-core exclusion between the tracer 
and a single obstacle. Note that for independent ob- 
stacles, only cr determines the void space allowing to 
consider the tracer as pointlike and the obstacles as hav- 
ing a radius cr. Then the obstacles naturally overlap re- 
sulting in clusters of excluded space. Alternatively, one 
may view the obstacles as pointlike and the tracer of ra- 
dius cr, reminiscent of a random Galton board. Figure [T] 
displays the case where the obstacles and the tracer are 
of the same size, with the highlighted particles moving 
through the frozen maze. 

A self-consistent mode-coupling kinetic theory has 
been introduced by Gotze, Leutheusser, and Yip [4011411 
l42l which provides a theory for various aspects of the 
dynamics of the Lorentz gas. In particular, their approx- 
imation gives an extremely accurate value for the first 
order density correction B3]| to the Lorentz-Boltzmann 
value for the diffusion coeflicient and correctly reflects 
the long-time tails in the velocity autocorrelation func- 
tion (VACF) Il44l . Furthermore, they correctly predict 
a localization transition at a critical obstacle density «* 
at which anomalous transport occurs. The competition 
between the critical law of the localization and the long- 
time tail then leads to density-dependent apparent expo- 
nents, a prediction that has been verified only recently 
by computer simulations fl3]| . Although the kinetic the- 
ory of the Lorentz model gives a remarkable description 
for low and moderate densities, the regime very close to 
the localization is more involved. 

The immediate vicinity of the arrest of transport is 
characterized by a diverging length scale describing the 
distribution of clusters of the void space. Already at 
moderate densities the void space falls apart into discon- 
nected finite components, hence confining the motion of 
the tracer particles. Thus even below the critical density, 
the system becomes non-ergodic in the sense that trac- 
ers initially placed in such a pocket remain there forever. 
Long-range transport occurs only on the infinite void 
space cluster percolating through the array of obstacles. 
At the transition point this infinite cluster ceases to exist 
and all tracer particles are trapped. Thus the origin of 
the localization is the underlying geometric percolation 
transition of the void space. Percolation itself has been 
the investigated for many years and a coherent picture 
has been established f32l, summarized below. 

Similar to a continuous phase transition, percolation 
exhibits a diverging length scale ^ ~ which cor- 
responds here to the linear dimension of the largest fi- 
nite pocket. Here the reduced obstacle density e = 
(n* - n*)/n* quantifies the distance to the critical den- 



sity and plays the role of a separation parameter The 
weight of the infinite void space cluster also vanishes as 
a power law, f „ ~ (-ef for e < 0, as the obstacle den- 
sity is increased. Directly at n*, the infinite void space 
becomes fractal with a fractal dimension c/f - d - /3/v, 
and in coexistence with the infinite cluster, there is a 
self-similar hierarchy of finite clusters. In contrast to 
thermal phase transitions a second length scale, viz. the 
mean-square cluster size {, plays an important role for 
the percolation transition. Using scaling arguments, one 
can show that its divergent behavior is determined by 
the same two exponents as ^ ~ \e\^'^'^^^^. The criti- 
cal exponents for percolation are known from computer 
simulations [32 , 38l and attain the approximate values 
V ^ 0.88 and /3 * 0.41 in three-dimensional space, 
c/ = 3. 

3. Motion on the infinite cluster 

Let us discuss first the dynamics of tracer particles 
that move on the percolating cluster only. Transport 
on the percolating cluster is expected to become frac- 
tal precisely at the critical obstacle density. After an 
elapsed time f a particle has moved typically a dis- 
tance f where (i„ is referred to as the walk dimen- 
sion and attains the value of d„ » 4.81 1381 in the 
three-dimensional Lorentz model. For the mean-square 
displacement 5r^(f) := ([R(f) - R(0)]^)^ this implies 
'5'"m(0 ~ f^^''", where the index oo indicates that aver- 
aging is only for particles initially located on the infi- 
nite cluster Below the critical obstacle density, trans- 
port should follow this subdiffusive behavior for inter- 
mediate times and then cross over to a diffusive motion 
6r-^ - 6Dcat for length scales larger than ^ ~ {-ey^, 
where the cluster appears homogeneous. Both aspects 
can be combined in the scaling law 

Jri(f;e) = f2/^»5ri(f), t^t/t, (1) 

i.e., all mean-square displacements have the same shape 
in a double-logarithmic representation. The dependence 
on the reduced obstacle density e is solely by a change 
of scale tx - t^ie). The scaling function 5r^(-) attains 
a constant for short rescaled times f <k 1, reflecting 
the anomalous transport. The length scale where the 
crossover to diffusion occurs is given by the correla- 
tion length, which then in turn determines the crossover 
time to tea ~ ^''^ ~ lei"''''"'. The crossover to diffu- 
sion is recovered by imposing dr^{t) ~ f'^^/rfw ^j- jQjjg 
rescaled times f » 1. As a consequence the corre- 
sponding diffusion coefficient vanishes as a power law, 
£>oo ~ f?'''""' ~ (_e)v(^fw-2)^ Since the particles on 
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Figure 2: Double-logarithmic representation of the mean-square dis- 
placements Srl^(t) for particles on the infinite cluster Obstacle density 
increases from top to bottom. The straight line indicates anomalous 
transport with a walk dimension of c/„ = 4.81. 



the infinite cluster are the only ones that contribute to 
long-range motion, one derives the scaling behavior of 
the diffusion coefficient for an all-cluster-average by 
Die) - Poaie)Doaie) ~ (-e)^ with the conductivity ex- 
ponent 11-13 + v{d^ - 2). 

We have performed extensive computer simulations 
for the Lorentz gas with ballistic particles as tracers. 
In this case the trajectory in configuration space con- 
sists of a series of straight lines connecting the scatter- 
ing events with the obstacles. We have employed de- 
terministic specular scattering which in particular con- 
serves the magnitude of the velocity v = |v(f)l- The tra- 
jectories are generated using the standard event-driven 
algorithm already used by Bruin [^461 combined with 
the usual method B7l l48l for calculating correlation 
functions online, optimized for exponentially large time 
scales. The infinite cluster has been identified by a stan- 
dard Voronoi tessellation to generate trajectories for the 
infinite-cluster-only averages. 

The simulation results for the mean-square displace- 
ment dr^(t) are displayed in Fig. |2] for various obsta- 
cle density below the localization transition. All curves 
show ballistic motion 6r^{t) - v^f for times smaller 
than the inverse collision rate or mean collision time 
Tc- In the low-density regime, n — > 0, the collision 
time grows, t^. ~ l/ncr^"'y, yet close to the percola- 
tion transition, the density dependence can be ignored. 
For low-obstacle densities the motion crosses over di- 
rectly to a linear increase of 6r^(t) in time, as is ex- 
pected for normal transport. As the critical density n* is 
approached a window of subdiffusive transport emerges 
and the mean-square displacements increase as a power 
law in accordance with theoretical prediction. Directly 
at the critical point n* ^ 0.839 the power law f^^''- is 



observed over more than 6 six decades in time. The 
long-time behavior for all densities below n* is again 
diffusive but sets in at later and later times as e t 0. A 
short inspection of the curves reveals that the shapes of 
the crossover functions are similar and scaling behavior 
is anticipated to hold on a semi-quantitative level. Yet 
including the corrections to scaling ll38l should corrob- 
orate the validity of the scaling hypothesis more gener- 
ally. A more detailed analysis of the scaling behavior of 
(5r^(f) can be found in P9l . 

4. All-cluster-averaged motion 

Conventionally the Lorentz model considers averages 
over all allowed starting positions of the tracer, i.e., 
both on any of the finite clusters as well as on the in- 
finite cluster. Long-range transport occurs only on the 
percolating cluster, but close to the threshold the finite 
clusters become arbitrarily large and contribute to trans- 
port at all scales. The motion of an ensemble of trac- 
ers in such a fractal landscape is also self-similar in the 
time-domain introducing a new exponent z characteris- 
tic of the subdiffusive motion at the critical point. The 
unconstrained mean-square displacement is defined as 
6r^{t) :- {[R{t) - R(0)]^^, where the average includes 
initial positions in the infinite and in finite clusters ac- 
cording to their weight as well as over the disorder. It 
increases again as a power law, 5r^(f) ~ f^^^, where 
the dynamic exponent z is larger than the walk dimen- 
sion c/w Since the all-cluster averaged mean-square 
displacement has akeady been discussed in Refs. [37] 
andHni here we focus on the time-dependent diffusion 
coefficient 



Z)(f):= ^<[R(0-R(0)]-v(0)) 

a 



(2) 



i.e., the correlation of the time-dependent displacement 
R(f) - R(0) with the initial velocity v(0) = R(0). The 
definition of D{t) holds in any dimension d, in particu- 
lar for d - 3. Because of time-translational invariance 
the time-dependent diffusion coefficient is related to the 
mean-square displacement via 



D(t) 



1 d 
2d dt 



(3) 



Similarly, the velocity autocorrelation function Z(f) :— 
{v(0 • v(0)) encodes the time-dependent motion 
via D{t) - Z(f') At' consistent with the Green-Kubo 
relation for the diffusion constant D :- D(t — > oo) as 
the long-time limit of the time-dependent diffusion co- 
efficient. 
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Figure 3: Double-logarithmic representation of the time-dependent 
diffusion coefficient Z)(0 = (l/6)d(5r^(f)/df. Density increases from 
top to bottom. The full line indicates subdiffusive behavior D{t) ~ 
p.lz-1 2[ ffig critical obstacle density x 0.839. The thick line in- 
dicates the power law r^l- characterizing the long-time decay in the 
localized phase. 



The simulation data for the time-dependent diffusion 
coefficient D{f) are displayed in Fig. [3] At short times, 
D{t) - v^t/d increases linearly and the curves for all 
scatterer densities superimpose. This reflects the free 
flight of the tracer until the first scattering event occurs 
at the collision time r^. At time scales of the order of the 
collision time, the time-dependent diffusion coefficient 
displays a maximum due to the cage eff'ect familiar from 
the zero crossing in the velocity autocorrelation func- 
tion. For longer times, D(t) decreases monotonically, 
effects of inertia become less and less relevant, and the 
structural relaxation dominates the transport. 

For low and moderate obstacle densities, D(t) rapidly 
saturates at its long-time limit D{t oo), which corre- 
sponds to the diffusion constant D of the tracer. Upon 
increasing the number of scatterers towards the per- 
colation threshold, diffusion is rapidly suppressed. In 
this regime D(t) displays a crossover from power-law 
decrease to saturation and one anticipates that scaling 
holds. A detailed analysis of the scaling properties of 
the three-dimensional Lorentz model in terms of the 
mean-square displacement can be found in Ref. i37| the 
scaling of the time-dependent diffusion coefficient for a 
two-dimensional Lorentz model is discussed in Ref.lSTI 

Precisely at the percolation threshold, D(t) ~ fi^^^^ 
is expected to hold for f — > oo as is inferred from the 
corresponding subdiffusive behavior of the mean-square 
displacement. Our simulations display such anomalous 
transport over more than five decades in time allow- 
ing the determination of the exponent z = 6.25, con- 



sistent with the prediction that the universality class of 
the three-dimensional Lorentz model coincides with the 
one of a random resistor network with power-law dis- 
tributed weak conductances Il371l50ll . 

In the high-density regime, n* > n*, the time- 
dependent diffusion coefficient vanishes in the long- 
time limit reflecting the fact that all tracers are trapped 
in finite pockets of the void space. Again, close to crit- 
icality the data follow the anomalous behavior until a 
crossover time and then decreases more rapidly. Yet 
the long-time behavior is not characterized by an ex- 
ponential decay since power law distributed exit rates 
from cul-de-sacs lead to long-time tails in the local- 
ized phase [34|. For the three-dimensional case, the ve- 
locity autocorrelation function (VACF) is predicted to 
decay as r^^^ which corresponds to r^^^ for the time- 
dependent diffusion coefficient as indicated in Fig. [3] To 
the best of our knowledge, this algebraic decay has not 
been observed before. 

5. Frequency-dependent conductivity 

In the context of ion conductors one is interested in 
the frequency-dependent complex conductivity cr(w) 
which quantifies the alternating electric current den- 
sity as response to an frequency-dependent ho- 
mogeneous electric field E(w) in the linear regime as 
j(cL>) - cr(a))E(tij). For disordered materials that are 
statistically isotropic the current is parallel to the elec- 
tric field and the conductivity o-{a)) transforms as a 
scalar. A decomposition into real and imaginary parts, 
cr(aj) - Re[cr(w)] + 1 Im[cr((y)], reveals that the in-phase 
component Re[cr(fa;)] > describes the loss due to fric- 
tion, whereas Im[cr((y)] encodes the storage of energy in 
the system. For the important case that the conducting 
ions can be viewed as independent, the conductivity is 
obtained by cr(w) = <7*nionju(w), where q* is the effective 
charge of the ions, nio„ their number density, and fJ.(cj) 
the frequency-dependent mobility. 

The linear response theorem extends the well-known 
Einstein relation D = fiksT to finite frequencies Z(a>) - 
kBTjiiiij), where the frequency-dependent generaliza- 
tion Z{<jj) of the diffusion coefficient D is given by 

1 r°° 

Z{(jS) = - <v(f) • v(0)) exp(iwf)df , (4) 
a Jo 

i.e., as the one-sided Fourier transform of the veloc- 
ity autocorrelation function. In particular, the Green- 
Kubo relation for the diffusion coefficient is recovered 
for the stationary case, w = 0. By partial integration one 
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Figure 4: Frequency-dependent diffusion coefficient Re[Z((i))] on 
double-logarithmic scales for various obstacle densities below tlie lo- 
calization threshold. Densities increase from top to bottom. The bro- 
ken line indicates the power law w'"-'- expected for anomalous trans- 
port at the critical density. 



Figure 5: Frequency-dependent susceptibility Re|;f(ai)] on double- 
logarithmic scales for various obstacle densities above the localiza- 
tion threshold. Densities increase from top to bottom. The brolcen 
line indicates the power-law divergence (j^^/: expected for anoma- 
lous transport at the critical density. 



obtains a representation of Z(w) in terms of the time- 
dependent diffusion coefficient 

r*oo 

Z(w) = D-iw I [D(f)-£>y"'df, (5) 
Jo 

with the diffusion coefficient D - D{t — > oo). 

Figure |4] displays the real part of the frequency- 
dependent diffusion coefficient, Re[Z((w)], obtained via 
a numerical one-side Fourier transform from our sim- 
ulation data according to Eq. |5] At the critical density 
n*, the power-law dependence Re[Z(w)] ~ w'"^^" holds 
over almost five decades in frequency. The fractal be- 
havior of Z{(L)) is inherited from an algebraic long-time 
decay of of the VACF, Z(f) ~ f^^"""^, or equivalently from 
the subdiffusive increase of the mean-square displace- 
ment, 6r-(t) ~ f-l^. 

For obstacle densities below n*, the curves attain a 
finite value at low-frequencies which is identified with 
the diffusion constant D - Zio) = 0). At high frequen- 
cies, the frequency-dependent diffusion coefficient ap- 
proaches the critical law. The two regimes of anoma- 
lous transport and diffusion merge at a characteristic 
crossover frequency u)^ which is shifted to lower values 
as the critical density is approached. 

These observations suggest the scaling behavior 

Z(<y; e) = (-\(jS)^'^^~X±{oS) , Gj - to/cj^- , (6) 

in the close vicinity of the localization transition. Here 
the scaling functions and Z-i-) refer to densities 

above {+) and below (-) the critical density, respec- 
tively. For large rescaled frequencies w » 1, the tracers 



cannot explore the entire landscape and transport is frac- 
tal, -3±((y — > oo) — > const, i.e., indistinguishable from 
the critical point. In particular, the constant is the same 
on both sides of the transition. 

On the conducting side, Z-icj) ~ (-lo;)^^'"' for cb — > 
is enforced by the requirement of a finite conductivity, 

1—2/" 

which in turn implies the scaling law D(e) ~ a>^- '~. 
Defining the conductivity exponent fi via D(e) ~ (-e)'' 
yields to, ~ (-e^'^'-^K 

On the insulating side, it is advantageous to discuss 
the complex polarizability = cr(w)/(-i(y) rather 

than the conductivity itself. Let us define the sus- 
ceptibility ^(w) = Z(a))/(-ia)) by dividing out factors 
that are irrelevant for the present discussion, x(^) = 
(q*nion/kBT)x{a>). Then^(a;) has the dimension of the 
square of a length. A finite polarizability is obtained by 
imposing 2^+(<y) ~ {-id>)^^~ for low frequencies w «: 1. 
The static susceptibility diverges upon approaching the 
critical obstacle density ^(w = 0, e) ~ . Be- 

cause of geometrical reasons this susceptibility should 
diverge in the same way as the mean-square cluster size 
Xioj ^ Q,e) ~ f ~ e-(2v-«. Combining the pre- 
ceding arguments reveals the scaling relation between 
the conductivity exponent and dynamic exponent z to 
Z = (2v - P + fj)/(v - /3/2). Then one also observes 
that the crossover frequency scales as ~ f^', i.e the 
inverse of the crossover time f , of the mean-square dis- 
placement for the infinite cluster only. 

The simulation results for the susceptibility on the 
insulating side are displayed in Fig. [5] As discussed 
above, the real part becomes finite in the static case 
reflecting the finite polarizability of the system. Ap- 
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Figure 6: Self-intermediate scattering function Fsiq, t) on a logarith- 
mic time scale at the critical density n* = 0.839. Wavenumbers in- 
crease from top to bottom. The full lines correspond to the simulation 
data, the broken ones to the Gaussian approximation F'^'^'^'^^iq, t), and 
the dotted ones include the first non-Gaussian correction F^'^^iq, t). 



preaching the delocahzation transition by diluting the 
obstacles gives rise to a strong increase of the static po- 
larizability. At the critical density, the susceptibility dis- 
plays power-law behavior in the frequency with finite 
response only for alternating external electric fields. 

6. Self-intermediate scattering function 

Spatio-temporal information on the dynamics of the 
tracer can be extracted by considering the self- 
intermediate scattering function (ISF) defined by 



F,(^,f):=<exp(iq.AR(f)> 



(7) 



Thus Fsiq, t) is the characteristic function of the dis- 
placements AR(f) := [R(f) - R(0)] considered as the 
random variable. We have anticipated already that due 
to statistical isotropy Fsiq, does not depend on the 
direction of the wavevector q, but only on its magni- 
tude q - |q|. The intermediate scattering function en- 
codes all moments of the displacement via derivatives, 
for instance the mean-square displacement can be ob- 
tained as 5r^(f) = -d\imq^od^Fsiq,t)/dq^. Scattering 
techniques such as neutron-spin echo or photon corre- 
lation spectroscopy have direct access to the ISF and 
there hq plays the role of the momentum transfer from 
the tracer to the neutrons or photons. The simulation re- 
sults for the ISF for small to moderate wavenumbers q 
for obstacle density directly at n* are displayed in Fig.|6] 
As expected transport becomes slower on larger length 
scales. In contrast to glassy dynamics characterized by 
a two-step relaxation process, here no plateau develops 



at intermediate times. Rather all intermediate scattering 
function appear to saturate at some finite level. Since 
the infinite cluster has non-extensive weight at the criti- 
cal point ||32]| its dynamics does not contribute to the ISF 
shown. Yet the infinite hierarchy of finite clusters coex- 
isting with the percolation cluster leads to non-ergodic 
behavior as manifested in the ISF, see below. Since the 
finite clusters are present even below the critical den- 
sity, one should in fact expect non-ergodic behavior due 
to the falling apart of the configuration space at all finite 
densities. 

The spatio-temporal information encoded in the ISF 
is revealed by gradually changing the wavenumber q. 
To highlight the non-trivial behavior we have included 
the Gaussian approximation id = 3) 



F?--(^,f) = exp(-^W(f)/6), 



(8) 



in Fig. |6] This approach ignores all cumulants except 
for the second, which reduces to the mean-square dis- 
placement Sp'it). The Gaussian model accounts approx- 
imately for the increase of time scale upon changing 
the wavenumber and also provides a quantitative de- 
scription of the initial decay of the ISF. Interestingly 
the Gaussian model becomes worse on increasing the 
length scale, demonstrating the fact that transport at 
the critical density never becomes Gaussian even at 
the largest scales. A similar observation on the planar 
Lorentz model has been made recently in the context of 
fluorescence correlation spectroscopy |,52J . Including 
also the first non-Gaussian parameter (NGP), 



aiit) 



[ARitff) I (AR(f)2) 



1, 



(9) 



in the intermediate scattering function. 



Ff°''iq,t)^C 



1 + -Q'2(0 



q^dr^t) 



(10) 

describes the simulation data even to times larger by a 
factor of 10. Nevertheless, this approximation too fails 
on a qualitative level for long times since it predicts a 
decay of all ISFs to zero due to the unbounded increase 
of the mean-square displacement. 



7. Lamb-Mossbauer factor 

The long-time limit of the intermediate scattering func- 
tion 

fsiq):^limFsiq,t), (11) 
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is referred to as Lamb-Mossbauer factor or non- 
ergodicity parameter, also known as the Edwards- 
Anderson parameter in the spin glass community [IJ. A 
non-vanishing fsiq) indicates that dynamic correlations 
are persistent forever implying that the dynamics is non- 
ergodic. Roughly speaking fAq) measures the fraction 
of particles that are trapped on length scale 2nlq. For 
the Lorentz model the presence of finite clusters at all 
densities suggests that the non-ergodicity transition has 
to be disentangled from the localization transition. A 
clear distinction can be made by considering the long- 
wavelength limit /j(0) = lim^^o/j(^) of the Lamb- 
Mossbauer factor, since 1 - /j(0) = P{e) where P(e) de- 
notes the weight of the percolating cluster. Since above 
the localization threshold n* > n* the infinite cluster 
ceases to exist, fs{Q) attains the value unity, whereas in- 
creasing the obstacle density from below leads to an ex- 
pected critical behavior 1 - /5(0) ~ {-£)^- 

For finite wavenumbers this sharp transition is 
rounded off" and crossover scaling should hold 1361 . A 
convenient form is obtained by postulating 1 - fsiq', e) = 
\£ff±(<lO, where the scaling functions f± refer to den- 
sities above the threshold (+) and the below the lo- 
calization transition (-). The critical behavior of the 
long-wavelength limit is recovered by requiring — > 
0) = const., and — > 0) = 0, respectively. For 
large wavenumbers and densities close to the transition 
the dynamics cannot resolve the finiteness of the cor- 
relation length and should become independent of the 
separation parameter, which is achieved by imposing 
f±ix) ~ x^^"". This in turn implies a singular depen- 
dence of the Lamb-Mossbauer factor on the wavenum- 
ber 1 - fsiq; e - 0) ~ q^l'' directly at the threshold. 

The Lamb-Mossbauer factors obtained from our sim- 
ulation results for the ISF as extrapolation to infinite 
times are displayed in Fig. [7] Upon increasing the den- 
sity fsiq) increases monotonically and the shape devel- 
ops into a rapid crossover for small wavenumbers. Yet 
for finite wavenumber nothing peculiar happens at the 
percolation threshold. Only the extrapolation to the long 
wavelength limit indicates the localization transition in 
the sense that it converges to unity rapidly as n* is ap- 
proached from below. Due to the large statistical un- 
certainties we have not attempted to validate the scaling 
law. Nevertheless the approach of the long-wavelength 
limit towards unity is compatible with the power law 
expected for the weight of the infinite cluster. 

8. Summary and Conclusion 

Transport in disordered material as studied in terms 
of the Lorentz model displays many facets that go be- 



n*. 



^ q = 2n/30cr 
^ q = 2nlWa- 
■ q — 27r/3cr 









0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 



Density n* 

Figure 7: Lamb-MoBbauer factor f,,{q) = Fs(q, ? — > oo) as a function 
of the obstacle density n*. As the wavenumber decreases the curves 
steepen. The extrapolation to zero wavenumber exhibits a sharp tran- 
sition at the critical obstacle density «* . The broken line corresponds 
to a power law {-ef and serves as a guide to the eye. 

yond the simple subdiffusive increase of the mean- 
square displacement. We have elucidated some of these 
aspects as they are manifest in the mean-square dis- 
placement for tracers that are confined to the percolating 
cluster only, the time-dependent diff'usion coefficient for 
an all-cluster average, the frequency-dependent conduc- 
tivity and susceptibility in response to an AC electric 
field, as well as the intermediate scattering function and 
its long-time limit, the Lamb-Mossbauer factor. 

The percolation transition of the underlying void 
space induces a self-similar geometric fractal with 
anomalous transport at the localization threshold. The 
interplay of time and length scales leads to a striking 
behavior of the complex transport properties many as- 
pects of which can be rationalized in terms of scaling 
laws. The scaling predictions require three indepen- 
dent exponents as input and we have demonstrated that 
all three values can be revealed by considering differ- 
ent limits of certain two-time correlation functions. In 
particular, the detailed study shows that simplified con- 
cepts such as fractional Brownian motion or continuous 
time random walks (CTRW) cannot explain the dynam- 
ics as soon as spatio-temporal information beyond the 
mean-square displacement at the threshold is studied. 

The localization transition in the Lorentz model in 
its simplest variant has now been studied from vari- 
ous aspects and nice agreement of extensive computer 
simulations with theoretical notions has been achieved. 
An open question is how anomalous transport occurs in 
a class of glass-forming liquids that displays a strong 
separation of time scales, and a possible route to an 
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answer is to include more microscopic details into the 
Lorentz model, thereby gradually adding more com- 
plexity. These changes could include studying the mo- 
tion of a tracer in a frozen matrix of correlated obsta- 
cles generated from a snapshot of a real liquid. In par- 
ticular, one has to clarify the role of the narrow gaps 
that emerge in the frozen labyrinth and their influence 
on the value of the dynamic or conductivity exponent. 
Furthermore, it is interesting to investigate the effect of 
interacting tracers inside the quasi-arrested host matrix 
since for long times the tracers will strongly interact in 
the ramified network of channels. A further interesting 
aspect is the interaction of a dynamic maze on the par- 
ticle, i.e., the effect of vibrations of the slow component 
on the dynamics of the tracer. 
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